Distributed delays stabilize neural feedback systems 
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We consider the effect of distributed delays in neural feedback systems. The avian optic tectum 
is reciprocally connected with the nucleus isthmi. Extracellular stimulation combined with intra- 
cellular recordings reveal a range of signal delays from 4 to 9 ms between isthmotectal elements. 
This observation together with prior mathematical analysis concerning the influence of a delay dis- 
tribution on system dynamics raises the question whether a broad delay distribution can impact the 
dynamics of neural feedback loops. For a system of reciprocally connected model neurons, we found 
that distributed delays enhance system stability in the following sense. With increased distribution 
of delays, the system converges faster to a fixed point and converges slower toward a limit cycle. 
Further, the introduction of distributed delays leads to an increased range of the average delay value 
for which the system's equilibrium point is stable. The enhancement of stability with increasing 
delay distribution is caused by the introduction of smaller delays rather than the distribution per se. 



PACS numbers: 87.19.L-, 87.18.Sn, 87.10.-c 

I. INTRODUCTION 

The signal flow in the brain is not just feedforward; 
rather, feedback dominates most neural pathways [1[ . Of- 
ten pairs of reciprocally connected neurons are spatially 
separate by several millimeters. For instance, the primate 
corticothalamic feedback loop extends over a distance of 
approximately 100 mm. Thus, for a typical action po- 
tential speed of 1 mm/ms we expect a signal delay of 
100 ms. When signal delays are larger than the neural 
response time, complex loop dynamics emerge 0,(3, H[. 

For reciprocally connected populations of neurons, 
large delays can introduce another dimension, namely the 
distribution of delay times. Such a distribution could be 
an epiphenomenon in the evolution of larger brains, or 
it could be of adaptive significance. Work from applied 
mathematics states an influence of the distribution of de- 
lay times on system dynamics @, H, H B B 03 ■ Intrigued 
by the latter possibility, we asked two questions: What is 
the distribution of delay times in an experimentally ac- 
cessible neural feedback system? What is the impact of 
distributed delays on a mathematically tractable neural 
model feedback system? 

We measured the distribution of delay times in the 
isthmotectal feedback system of birds [Fig. QJa)] (Iil.[l2l]. 
The avian isthmic nuclei (parabigcminal nucleus in mam- 
mals) receive a topographically organized projection from 
the tectum (superior colliculus in mammals), to which 
they project back and have been conjectured to medi- 
ate spatiotemporal attentional mechanisms [HI, [lj, [l5| . 
The isthmic nuclei in birds consist of three substructures: 
pars parvocellularis (Ipc), pars magnocellularis (Imc), 
and pars semilunaris (SLu) that are spatially separated 



'Electronic address: rw@physics.wustl.edu 



from the tectum [lfj, [17] . In response to visual stimu- 
lation, the Ipc neurons undergo a transition from quies- 
cence to rhythmic firing [Tol . [lg | . Delays can drive a neu- 
ral feedback loop over a stability boundary resulting in 
oscillatory behavior [H, [20|, HI IH [H Hi] . To elucidate 
the impact of a delay distribution on the system dynam- 
ics, we investigated, through numerical simulations and 
mathematical analysis, a model of reciprocally coupled 
neurons with distributed delays. 



II. MEASURED DISTRIBUTION OF DELAYS 

To measure the signal delays between pairs of isth- 
motectal elements, we obtained intracellular whole-cell 
recordings from identified neurons in a midbrain slice 
preparation and stimulated groups of presynaptic neu- 
rons or axons with brief electrical pulses delivered extra- 
cellularly [Fig. [ljb)]. Neurons were identified by their 
location within the midbrain slice preparation and for a 
subset of recorded neurons we obtained additional iden- 
tification via intracellular fills [HI [l7| • 

A subpopulation of tectal layer 10 (L10) neurons 
project to both the ipsilateral Ipc and Imc in a topo- 
graphic fashion 0, [l7|, [H [H, US HI- Their apical 
dendrite courses straight up to layer 2 with few ramifi- 
cations, and basal dendrites reach down to the border of 
layer 13. Retinal axon terminals overlap with the api- 
cal dendrite in tectal layers 2 to 7 [29|, y(| ■ We placed 
a stimulus electrode in layer 2 to 4 (L2-4) and recorded 
from L10 neurons with whole-cell recordings in response 
to L2-4 stimulation. The delays from the beginning of the 
stimulus pulse to the onset of the L10 response ranged 
from 4 to 15 ms with a mean delay of 6.9 ms and a stan- 
dard deviation (SD) of 1.3 ms (n = 15 cells) [Fig. Ufa)]. 
Tectal L10 neurons are a heterogeneous population p7j . 
Therefore, only filled L10 neurons with axons originating 
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FIG. 1: (Color) Schematic of the isthmotectal circuitry 
and representative response to electrical stimulation, (a) 
Schematic of the isthmotectal circuitry. RGC axons (black 
arrows) enter in upper tectal layers. A subpopulation of tectal 
L10 neurons (red) projects to Imc and Ipc. The Imc nucleus 
consists of two populations of neurons (blue); one project- 
ing broadly back to lower tectal layers and one projecting 
broadly to the Ipc nucleus. Ipc neurons (green) project back 
to the tectum with axons reaching into upper tectal layers, 
(b) Intracellular recording from an Ipc neuron in response to 
electrical stimulation in tectal L2-4. 



from the dendrite were included in this analysis. Since 
L10 neuron dendrites can reach up to L2, the possibility 
of unwanted direct electrical, rather than synaptic, stim- 
ulation of L10 neuron dendritic endings arises. At the 
end of a recording session, we evaluated the nature of 
stimulation by blocking chemical synaptic transmission 
via the block of Ca-channels by replacing Ca 2+ in the 
saline with Mg 2+ . 

We measured signal delays between optic tectum and 
individual Ipc neurons via RGC axon stimulation or L10 
neuron dendrite stimulation, with a stimulus electrode 
placed in tectal L2-4. In the first case, the group of stimu- 
lated RGC axons stimulates a population of L10 neurons, 
which in turn stimulates a large number of Ipc neurons. 
In the second case, L10 neurons are stimulated directly. 
This stimulus paradigm provided a high chance of record- 
ing from an Ipc neuron that received tectal synaptic in- 
puts. The delays from the beginning of the stimulus pulse 
to the onset of the Ipc neuron response ranged from 5 to 
19 ms (n — 17 cells) [Fig. [H[b)]. As expected from the 
stimulus paradigm, the distribution of delays is bimodal. 
We suspect that the first bump (5 to 9 ms range) is domi- 
nated by direct L10 dendrite stimulation (mono-synaptic 



pathway LlO-Ipc); whereas the second bump (11 to 19 ms 
range) is dominated by RGC axon stimulation, which ini- 
tiates the bi-synaptic pathway RGC-LlO-Ipc. From the 
first bump in the histogram we estimate a mean delay of 
6.5 ms and a SD of 1.4 ms for the mono-synaptic pat hway 
LlO-Ipc. Since Ipc neuron axons can reach up to L2 [17j |. 
the possibility of unwanted direct electrical stimulation 
of Ipc axons arises. At the end of a recording session, we 
evaluated the nature of stimulation by blocking chemical 
synaptic transmission via replacing Ca 2+ in the saline 
with Mg 2+ . 

Using a stimulus paradigm similar to the one described 
above, we measured signal delays between L10 and in- 
dividual Imc neurons. We placed a stimulus electrode 
in L2-4 for stimulation of RGC axons or L10 neuron 
dendrites and recorded from Imc neurons with whole- 
cell recordings in response to L2-4 stimulation. The sig- 
nal delays ranged from 4 to 19 ms (n — 17 cells) and 
the distribution was bimodal [Fig. HJc)]. As described 
above, the first bump is likely to be dominated by the 
mono-synaptic pathway (L 10- Imc), whereas the second 
bump is likely to be dominated by the bi-synaptic path- 
way (RGC-LlO-Imc). The first bump in the histogram 
yielded a mean delay of 5.2 ms and a SD of 0.9 ms. Since 
Imc axons terminate in tectal layers 10 to 13 [l||, the 
possibility of direct Imc axon stimulation via stimulus 
electrodes in L2-4 does not arise. 

The Imc nucleus consists of two cell types, one of which 
projects to the Ipc nucleus with a broad and dense pro- 
jection of axonal arbors [lfl [HI, [HJ. We positioned a 
stimulus electrode in the Imc nucleus and recorded from 
Ipc neurons with whole-cell recordings in response to Imc 
stimulation. The signal delays ranged from 3 to 8 ms with 
a mean delay of 4.3 ms and a SD of 1.1 ms (n = 12 cells) 
[Fig- Eld)]. Care had to be taken about the interpreta- 
tion of the Imc stimulation experiments. The stimulus 
electrode in the Imc nucleus stimulates 4 elements: L10 
neuron axons, Ipc neuron axons passing through the Imc 
nucleus, and two populations of Imc neurons; one pro- 
jecting to tectum and the other projecting to Ipc. To 
filter out the Imc to Ipc synaptic connection, we stim- 
ulated in an area of the Imc nucleus that did not cor- 
respond to the topographic location of the recorded Ipc 
neuron, thus avoiding both antidromic stimulation of the 
axon from the recorded Ipc neuron as well as avoiding or- 
thodromic stimulation of the L10 axons passing through 
the Imc nucleus on their way to the same location in 
the Ipc nucleus. At the end of a recording session, we 
applied bicuculline to verify that the synaptic inputs to 
the recorded Ipc neuron were indeed from the stimulated 
GABAergic Imc neurons (GABA: gamma-aminobutyric 
acid). The responses disappeared when 100 fiM bicu- 
culline was added to the bath (data not shown) thus 

(i) indicating that the responses were of synaptic origin 
(rather than antidromic Ipc or L10 axon stimulation) and 

(ii) confirming that GABA is the transmitter as had been 
suggested by anatomical studies 

The Ipc nucleus has topographical reciprocal connec- 



3 




Delay (ms) Delay (ms) Delay (ms) 

FIG. 2: (Color) Measured distribution of signal delays between isthmotectal elements and plot of the corresponding gamma 
distribution [red curves in (a) through (d)] with the same mean and standard deviation, (a) L2-4 to L10. (b) L2-4 to Ipc. 
(c) L2-4 to Imc. (d) Imc to Ipc. (e) Ipc to L10. (f) Imc to L10. 



tions with the tectum 0, EE EE EE El- The effer " 
ents from Ipc have large calibre axons and terminate in a 
columnar manner ranging from layers 2 to 12 [Fig. QJa)] 
03, EE EE EE E3, El- We applied local extracellular 
electrical stimulation of a group of Ipc neurons with a 
stimulus electrode placed in the Ipc nucleus. Such ex- 
tracellular electrical stimulation also stimulates L10 ax- 
ons antidromically. The fast L10 neuron antidromic re- 
sponses were distinguishable from the much slower and 
long-lasting synaptic responses. The additional direct 
activation of Imc axons in the Ipc nucleus does not in- 
terfere with this experiment, since the population of Imc 
neurons projecting to the Ipc nucleus is different from the 
population of Imc neurons projecting to the tectum. The 
yield for finding Ipc to L10 synaptic responses turned out 
to be very low. For the few cases we found, the delays 
ranged from 6 to 8 ms (n = 5 cells) [Fig. [2Je)]. 

The projection from individual Imc neurons to tectal 
layers 10 to 13 is broad and sparse [ljj]. We positioned a 
stimulus electrode in the Imc nucleus and recorded from 
L10 neurons with whole-cell recordings in response to 
Imc stimulation. The yield for finding Imc to L10 synap- 
tic responses turned out to be very low. For the two 
connected pairs we found, the signal delays were 3 and 
6 ms (n = 2 cells) [Fig. 12(f)]. The low yield and the 
interpretation of these experiments require some expla- 
nation. As mentioned above, a stimulus electrode in the 
Imc nucleus will stimulate four elements. To filter out 
the Imc to L10 synaptic connection, we stimulated in an 
area of the Imc nucleus that did not correspond to the 
topographic location of the recorded L10 neuron, thus 
avoiding both antidromic stimulation of the axon from 
the recorded L10 neuron as well as avoiding orthodromic 
stimulation of the Ipc axons passing through the Imc nu- 
cleus on their way to the same location of the tectum. 



At the end of a recording session, we applied bicucullinc 
to verify that the synaptic inputs to the recorded L10 
neuron were indeed from the stimulated GABAergic Imc 
neurons. For the two neurons, the responses disappeared 
when 100 fiM bicuculline was added to the bath (data not 
shown) thus indicating that the responses were of synap- 
tic origin; rather than antidromic L10 or orthodromic Ipc 
axon stimulation. 

In summary, these data show that the signal delays be- 
tween isthmotectal elements are distributed ranging from 
4 to 9 ms. 



III. DISTRIBUTED DELAYS AND THE 
DYNAMICS OF NEURAL FEEDBACK SYSTEMS 

What is the impact of distributed delays on a mathe- 
matically tractable neural model feedback system? To in- 
terpret the potential impact of the measured distribution 
of delays on the dynamics of neural feedback systems, 
we investigated a model system of two coupled Hopfield 
neurons 0, E2, EE Ell > described by the first-order delay 
differential equations 

duiit) 

— — = -Ui(t) +aitann|ti2(t — T2)\ , 

- -tta(t) + aatanh[«i(t-n)]. (1) 

at 

Here, u\(t) and U2(t) denote the voltages of the model 
neurons and t\ and t-i are the temporal delays, while a\ 
and 0,2 describe the coupling strength between the two 
neurons. In this analysis, the time variable is dimension- 
less. Translation to real time can be made by multiplying 
the dimensionless time variable with a membrane time 
constant, RC. The system of delay differential equations 
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FIG. 3: (Color) Mean delays and attractors. (a), (b) Dynam- 
ics of the two-neuron model system for gamma distributions 
with mean delay values of T = 0.7 [(a), fixed point] and [(b), 
limit cycle], respectively. For both cases, the standard devi- 
ation is 0% (blue), 25% (black), and 50% (red) of the mean 
delay. The initial condition is «i(i) = 0.30 and U2(t) — —0.28 
for — r < t < 0. (c) Gamma distribution for a mean de- 
lay value of T = 0.7 and a standard deviation of 0% (blue), 
25% (black), and 50% (red) of the mean delay, (d) Critical 
mean delay, To, where the Hopf bifurcation takes place, plot- 
ted against variance, (e), (f) Time constant for reaching the 
attractor for T = 0.7 (fixed point) and T = 2.0 (limit cy- 
cle), respectively, plotted against the variance of the gamma 
distribution. 



has a trivial stationary point at the origin, u\ — u 2 — 
[Fig. [3Ja)]. For a\a 2 < —1, the fixed point at the ori- 
gin is asymptotically stable as long as the mean of the 
time delays (ti + t 2 )/2 does not exceed a critical value 
r IH, If: 



n + t 2 
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The critical value To is determined by combinations of the 
product of the couplings alone [Eq. ©]. For couplings 
of opposite signs (e.g., a\a 2 < —1) and when the delays 
are increased, the origin becomes unstable and a limit 
cycle emerges via a supercritical Hopf bifurcation at (t± + 
t 2 )/2 = t [Fig. [3jb)]. The critical value, To , decreases 
with decreasing value of the product of the couplings a\a 2 
below —1. In other words, oscillations can be achieved by 
either increasing the delays or by increasing the absolute 
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FIG. 4: (Color) Dynamics of the two-neuron model system 
with discrete delays, (a) Dynamics of the system with two 
fixed delays of 0.1 and 0.7 (green); one fixed delay of 0.1 
(red); and one fixed delay of 0.7 (black), (b) Distance D(t) = 
\J u\(t) + u? 2 (t) from (0, 0) vs. time for the above cases. 



value of the coupling strengths of opposite signs. 

For a distribution of delays we replace the coupling 
term in JT]) with a weighted sum over similar terms but 
with different delays 

du\(t) 

dt 
du 2 (t) 

dt 

The delay kernel £(r) is normalized to satisfy 
g?t£(t) = 1 . For simplicity, we chose the delay ker- 
nels to be identical for both legs of the loop. We chose 
the delay kernel to be a gamma distribution, 



-ui(t) + ai / dr^r) tanh[u 2 (t — t)] , 
Jo 

-u 2 {t)+a 2 / drf(T)tanh[ui(t-T)].(3) 
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where T is the mean delay, v is the variance of the 
gamma distribution, and the gamma function is defined 



as rOr)= J?? 



dt. The gamma distribution was 



chosen because it has the biologically plausible feature 
to vanish for delays approaching [Fig. EJc)]. For the 
coupling strengths we chose a\ = — 2 and a 2 = 1 for all 
simulations. 

The parameters to vary are the mean delay, T, and the 
variance, u, of the gamma distribution. As these param- 
eters are changed, the fixed point at the origin changes 
from a stable fixed point to an unstable fixed point sur- 
rounded by a stable limit cycle and vice-versa (Hopf bi- 
furcation). This transition takes place when the roots, A, 
of the characteristic equation for the system ([3]) , 



A + 2A + 1 — a x a 2 1 



0. 



(5) 



are purely imaginary. The characteristic equation is ob- 
tained by demanding that the solution to ([3]) behaves as 
U\(t) — Ae xt , u 2 (t) — Be xt near the fixed point. Substi- 
tuting A = iu>, where tu is real, we have 



, , ILOV 

— lu + 2iu> + 1 — a\a 2 \ 1 H — — 



-2T 2 /v 



= 0. 



(6) 
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Separating real and imaginary parts, we get two equa- 
tions from which we can numerically eliminate to. How- 
ever, there are multiple solutions for this. For a given 
variance v, the solution with the minimum positive mean 
delay T, determines the critical mean delay To at which 
the fixed point at the origin loses its stability and a stable 
limit cycle emerges. Our analysis shows that the intro- 
duction of distributed delays (increasing variance) leads 
to a smaller limit cycle [Fig. [3jb)]. Furthermore, the 
critical mean delay T increases with increasing variance 
[Fig.Hd)]. 

To estimate the time constant for reaching an attrac- 
tor, we calculated the distance Dg(t) — y/u\{t) + 
from the origin along a given polar angle 9 in the u\-ui 
space. Assuming an exponential dependence, a fit of an 
exponential function to the simulated Dg(t) values pro- 
vided the time constant for that polar angle. We repeated 
the procedure for 360 polar angles in 1-deg increments 
and took the final time constant to be the mean of the 
360 time constants at given polar angles. This analysis 
shows that increasing variance makes the convergence to 
the fixed points faster [Fig. Efe)] and the convergence to 
limit cycles slower [Fig. E^f)]. 

In summary, distributed delays increase the parame- 
ter region with fixed-point behavior and accelerate the 
convergence to the fixed point. 

The enhanced fixed-point stability with distributed de- 
lays actually comes from the contribution of the smaller 
delays, rather than the distribution per se. To illustrate 
this insight, we compare the dynamics of three systems 
with delays given by one or two superimposed delta- 
distributions with peaks at: (i) r = 0.7, (ii) both r = 0.1 
and r = 0.7, and (iii) r = 0.1 (Fig.g|). The distributed 
system (ii) converges faster than system (i), but the dis- 
tributed system (ii) converges slower than system (iii). 
In other words, adding a longer delay of r = 0.7 to the 
t = 0.1 system slows the convergence, whereas adding a 
shorter delay of r = 0.1 to the r = 0.7 system acceler- 
ates the convergence. Thus, it is not the distribution of 
delays per se, but the contribution of shorter delays in 
the distribution that enhances fixed-point stability. 



IV. DISCUSSION 

For large brains with finite signal propagation veloci- 
ties, delays are a fact of life. In some feedforward path- 
ways, such as the vertebrate optic nerve, delays can be 
specific to the retinal ganglion cell type thus leading to 
differences of arrival time for different retinal represen- 
tations of the visual stimulus [13, [H, [H, H(| ■ In other 
feedforward sensory pathways, such as the avian nucleus 
laminaris, delays are used explicitly to evaluate interau- 
ral time differences [4l|. Delays in feedback loops play a 
fundamentally different role, as they can determine the 
dynamical behavior of the system [Z^, IH| . Specifically, 
for delays smaller than a critical value a neural feedback 
system may converge toward a steady-state, whereas for 



delays larg er than the critical value the system may os- 
cillate jUH^]. In nonlinear systems, the distribution of 
a system parameter can have unexpected effects on the 
systems dynamics 0, [47], EU . Consequently, if delay is a 
relevant parameter in neural feedback systems, as stated 
above, it is important to investigate the impact of delay 
distributions on the system dynamics. 

Parameters in biological system are usually distributed 
over some range. Therefore, these systems must be ro- 
bust in the sense that the parameter variability should 
not detain the biological system from functioning cor- 
rectly. Moreover, a system architecture in which the pa- 
rameter variability actually enhances the system's per- 
formance would be particularly desirable. In this study, 
we have quantified the distribution of delays in the avian 
isthmotectal feedback loop. Furthermore, by investigat- 
ing a mathematical model of coupled neurons with dis- 
tributed delays, we have demonstrated that distributed 
delays enhance the stability of the system, where the sta- 
bilizing effect arises from the contribution of smaller de- 
lays introduced through the delay distribution. Since the 
functional role of the isthmotectal feedback loop remains 
mostly unclear to date, it is not obvious whether this 
stabilizing effect is beneficial to the system's functioning. 
Further understanding of the neuronal processes in the 
isthmotectal feedback loop will be necessary to answer 
this question. 



V. DISCUSSION 

For large brains with finite signal propagation veloci- 
ties, delays are a fact of life. In some feedforward path- 
ways, such as the vertebrate optic nerve, delays can be 
specific to the retinal ganglion cell type thus leading to 
differences of arrival time for different retinal represen- 
tations of the visual stimulus [13, HE [H, S3] • In other 
feedforward sensory pathways, such as the avian nucleus 
laminaris, delays are used explicitly to evaluate interau- 
ral time differences [4l|. Delays in feedback loops play a 
fundamentally different role, as they can determine the 
dynamical behavior of the system [22, EI| . Specifically, 
for delays smaller than a critical value a neural feedback 
system may converge toward a steady-state, whereas for 
delays larg er than the critical value the system may os- 
cillate [44ri45l|. In nonlinear systems, the distribution of 
a system parameter can have unexpected effects on the 
systems dynamics [H, H3, HU . Consequently, if delay is a 
relevant parameter in neural feedback systems, as stated 
above, it is important to investigate the impact of delay 
distributions on the system dynamics. 

Parameters in biological system are usually distributed 
over some range. Therefore, these systems must be ro- 
bust in the sense that the parameter variability should 
not detain the biological system from functioning cor- 
rectly. Moreover, a system architecture in which the pa- 
rameter variability actually enhances the system's per- 
formance would be particularly desirable. In this study, 
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we have quantified the distribution of delays in the avian 
isthmotectal feedback loop. Furthermore, by investigat- 
ing a mathematical model of coupled neurons with dis- 
tributed delays, we have demonstrated that distributed 
delays enhance the stability of the system, where the sta- 
bilizing effect arises from the contribution of smaller de- 
lays introduced through the delay distribution. Since the 
functional role of the isthmotectal feedback loop remains 
mostly unclear to date, it is not obvious whether this 
stabilizing effect is beneficial to the system's functioning. 
Further understanding of the neuronal processes in the 
isthmotectal feedback loop will be necessary to answer 
this question. 



VI. EXPERIMENTAL METHODS 

White Leghorn chick hatchlings (Gallus gallus) of less 
than 3 days of age were used in this study. All pro- 
cedures used in this study were approved by the local 
authorities and conform to the guidelines of the National 
Institutes of Health on the Care and Use of Laboratory 
Animals. Animals were injected with ketamine (40 mg 
per kg, i.m.). Brain slices of the midbrain were pre- 
pared following published protocols [SiJ HH, [H, [53[ . 
Briefly, preparations were done in °C, oxygenated, and 
sucrose-substituted saline (240mM sucrose, 3 mM KC1, 
5 mM MgCl 2 , 0.5 mM CaCl 2 , 1.2 mM NaH 2 P0 4 , 23 
mM NaHC0 3 , and 11 mM D-glucose). After decapita- 
tion, the brains were removed from the skull, and the 
forebrain, cerebellum, and medulla oblongata were dis- 
carded. A midsagittal cut was used to separate the tectal 
hemispheres. The tectal hemispheres were sectioned at 
500 /xm on a tissue sheer (Vibroslice, Camden and VF- 
200, Precisionary Instruments) in either the transverse 
or the horizontal plane. Slices were collected in oxy- 
genated saline (120 mM NaCl, 3 mM KC1, 1 mM MgCl 2 , 
2 mM CaCl 2 , 1.2 mM NaH 2 P0 4 , 23 mM NaHC0 3 , and 
11 mM D-glucose) and kept submerged in a chamber 
that was bubbled continuously with carbogen (95% oxy- 
gen, 5% C0 2 ) at room temperature. The slice was then 
transferred to a recording chamber (RC-26G, Warner In- 
struments) mounted on a fixed stage upright microscope 
equipped with DIC optics (BX-51WI, Olympus). The 
slice was held gently to the bottom of the chamber with 
an anchor of nylon threads, and the chamber was per- 
fused continuously with oxygenated saline at room tem- 
perature. The cells in L10, Imc, and Ipc are visible with 
DIC optics. 

Local electrostimulation was achieved by inserting 
bipolar tungsten electrodes under visual control into ei- 
ther the upper tectal retinorecipient layers (2 to 4), layer 
5b, or the isthmic nuclei Ipc or Imc with a three-axis 
micromanipulator (U-31CF, Narishige). Electrodes were 
custom-built from 50-/xm diameter, insulated tungsten 
wires (California Fine Wire) that were glued together 
with cyanoacrylate and mounted in glass microcapillar- 
ies for stabilization. The wires protruded several hundred 



m from the capillaries, and the tips were cut at an angle. 
Stimulus isolators (Isolated Pulse Stimulator 2100, AM 
Systems) generated biphasic current pulses (20 - 200 /xA, 
500 fjs). 

Whole-cell recordings were obtained with glass mi- 
cropipettes pulled from borosilicate glass (1.5 mm OD, 
0.86 mm ID, AM Systems) on a horizontal puller (P- 
97, Sutter Instruments and DMZ Universal Puller, Zeitz 
Instruments) and were filled with a solution containing 
100 mM K-Gluconate, 40 mM KC1, 10 mM HEPES, 0.1 
mM CaCl 2 , 2 mM MgCl 2 , 1.1 mM EGTA, 2 mM Mg- 
ATP, pH adjusted to 7.2 with KOH. Electrodes were 
advanced through the tissue under visual guidance with 
a motorized micromanipulator (MP-285, Sutter Instru- 
ments) while constant positive pressure was applied and 
the electrode resistance was monitored by short current 
pulses. Once the electrode had attached to a membrane 
and formed a seal, access to the cytosol was achieved by 
brief suction. Whole-cell recordings were performed with 
the amplifier (Axoclamp 2B, Axon Instruments and SEC- 
05L, npi-electronic) in the bridge mode (current clamp). 
The series resistance was estimated by toggling between 
the bridge and the DCC (discontinuous current clamp) 
mode. The series resistance was compensated with the 
bridge balance. Analog data were low-pass filtered (4- 
pole Butterworth) at 1 kHz, digitized at 5 kHz, stored, 
and analyzed on a PC equipped with an PCI-MIO-16E-4 
and Lab View software (both National Instruments). 

Labeling of a subset of recorded neurons was carried 
out as described previously [5(1 HH, [H, IHj]. In brief, 
whole-cell patch recordings were obtained as described 
above. Additionally, the electrode solution contained 
0.5% Biocytin (w/v) to label the recorded neurons. Indi- 
vidual cells were filled intracellularly with 2 nA of posi- 
tive current over 3 minutes. After recording and labeling, 
slices were kept in oxygenated ACSF for an additional 
30 minutes and subsequently fixed by immersion in 4% 
paraformaldehyde in PB for at least 4 hours. Slices were 
then washed in phosphate buffer (PB, 0.1 M, pH 7.4) 
for at least 4 hours, immersed in 15% sucrose in PB for 
at least 4 hours and then immersed in 30% sucrose in 
PB for 12 hours, and resectioned at 60 /xm on a freezing 
microtome. The sections were collected in PB and the 
endogenous peroxidase blocked by a 15-minute immer- 
sion in 0.6% hydrogen peroxide in methanol. The tissue 
was washed several times in PB, and then incubated in 
the avidin-biotin complex solution (ABC Elite kit, Vector 
Labs) and the reaction product visualized with a heavy- 
metal intensified DAB protocol. Following several washes 
in PB, the 60 /xm-thick sections were mounted on gelatin- 
coated slides, dried, dehydrated, and coverslipped. Sec- 
tions were inspected for labeled neurons, and only data 
from cells that could unequivocally be classified accord- 
ing to published criteria [HI [TtJ were taken for further 
analysis. Cells were reconstructed at medium magnifica- 
tion (lOx to 20x) with a camera lucida on a Leica micro- 
scope and projected onto the 2D plane. 
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